Imaging the symmetry breaking of molecular orbitals in carbon nanotubes 
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Carbon nanotubes have attracted considerable interest for their unique electronic properties. They 
are fascinating candidates for fundamental studies of one dimensional materials as well as for future 
molecular electronics applications. The molecular orbitals of nanotubes are of particular importance 
as they govern the transport properties and the chemical reactivity of the system. Here we show 
for the first time a complete experimental investigation of molecular orbitals of single wall carbon 
nanotubes using atomically resolved scanning tunneling spectroscopy. Local conductance measure- 
ments show spectacular carbon-carbon bond asymmetry at the Van Hove singularities for both 
semiconducting and metallic tubes, demonstrating the symmetry breaking of molecular orbitals in 
nanotubes. Whatever the tube, only two types of complementary orbitals are alternatively observed. 
An analytical tight-binding model describing the interference patterns of n orbitals confirmed by ab 
initio calculations, perfectly reproduces the experimental results. 

PACS numbers: 78.67.Ch,73.22.Dj,68.37Ef 



Single-walled carbon nanotubes (SWNTs) are fasci- 
nating candidates for fundamental investigations of elec- 
tron transport in one-dimensional systems, as well as for 
molecular electronics [H, 0, Q . The control of their elec- 
tronic bandstructure using doping or functionalization is 
a major challenge for future applications. In this con- 
text, the molecular orbitals of SWNTs play a major role 
since they are fully involved in the transport properties 
and the chemical reactivity of SWNTs [J]. In particu- 
lar their spatial electronic distribution is a key factor in 
chemical reactions [H| and in adsorption process Q. In 
transport experiments, it is an important parameter for 
the contact quality between a tube and an electrode 0. 
A complete understanding and control of the molecular 
orbitals of nanotubes is therefore crucial for the develop- 
ment of future nanotubes based applications. 

Scanning tunneling microscopy and spectroscopy 
(STM/STS) are unique tools to measure local electronic 
properties of SWNTs and correlate them with their 
atomic structure. STS studies have confirmed the pre- 
dictions of the simple zone folding tight-binding model 
which relates the metallic or semiconductor character to 
the chirality of SWNTs @, @, In addition, an elec- 
tronic "pseudogap" on metallic tubes has been measured 
as predicted theoretically when curvature effects or inter- 
tube interactions are fully considered ll|. Interference 



patterns have been observed in STM images close to de- 
fects such as chemical impurities, cap ends or in finite 
length nanotubes [H, EH, EiL EH- These observations 
have been explained qualitatively from interference ef- 
fects between the Bloch waves [la, LL2l ■ However, the 
molecular orbitals corresponding to the Bloch states as- 
sociated to the Van Hove singularities (VHS) have been 



less investigated. It has been predicted that, for defect- 
free semiconducting tubes, the electronic states would 
display broken symmetry effects in agreement with STM 
topographic images H, 3, T3|. The specific signatures 
of such Bloch states in metallic and semiconducting tubes 
can only be evidenced using local differential conductance 
spectroscopy, but experimental data are still scarce. 

Here, we report on a systematic study of the local den- 
sities of states (LDOS) in defect-free SWNTs at various 
VHS using a low temperature STM. For semiconduct- 
ing tubes, direct images of the two first highest occupied 
molecular-orbitals (HOMO-1, HOMO) and lowest un- 
occupied molecular-orbitals (LUMO, LUMO+l) states 
have been obtained. A spectacular C-C bond asymmetry 
is observed, revealing the complementarity in the symme- 
try of the respective wavef unctions. For metallic tubes, 
we observe the splittings of the first VHS and the oc- 
currence of a pseudogap. As a consequence, the wave- 
function symmetry breaking remains, which is demon- 
strated by conductance images showing strong symme- 
try variations within a few meV. Symmetry analysis of 
the wavefunctions within an analytical tight-binding ap- 
proach and calculations of STM images using density- 
functional based calculations in the Tersoff-Hamann for- 
malism give a complete theoretical description of these 
findings. 

Semiconducting tubes We present first a detailed 
study of a semiconducting tube. From the topographic 
image (Figflh.), we estimate a chiral angle about 15±1° 
and a diameter of 1.4±0.1 nm, leading to (15,5) chiral 
indices for the studied nanotube. The dl/dV spectra 
measured on the three different C-C bond orientations 
(FigHb) exhibit four peaks at -0.78, -0.37, +0.37, and 




FIG. 1: (a) large scale constant current topographic image of a (15, 5) semiconducting tube, (b) topographic image at atomic 
scale recorded during the CITS measurement with labelling of the three different bonds (Vtias = 1 V, I = 0.2 nA). The scheme 
indicates notations used in the tight-binding discussion. f a and 9 a are the three C-C bonds vectors and their angles with the 
tube axis respectively, (c) STS measured at three different C-C bonds, (d)-(g): Energy-dependent conductance images of the 
tube at -0.78 , -0.37 , 0.37, and 0.78 V respectively. The arrow indicates the tube axis. 



+0.78 V, corresponding to four VHS denoted VHS-2, 
VHS-1, VHS+1 and VHS+2, respectively (FigHfc). The 
wavefunctions at the corresponding energies will be de- 
noted HOMO-1, HOMO, LUMO and LUMO+1 respec- 
tively. The most striking feature of the spectra is that the 
relative intensity of the VHS peaks is position-dependent. 
Indeed, comparing the spectra measured above bond 1 
and bond 2, the four peaks vary in an alternate way, 
i.e. VHS-2 decreases while VHS-1 increases, VHS+1 de- 
creases and VHS+2 increases. A similar variation is also 
observed between bond 2 and bond 3 but with a smaller 
amplitude. 

FigQJl-g show energy-dependent differential conduc- 
tance images (dl/dV maps) of the semiconducting tube. 
The significant deformation of the hexagonal lattice is 
mostly due to the drift occuring during the long acqui- 
sition time needed for current imaging tunneling spec- 
troscopy (CITS) measurements as well as distorsion in- 
duced by the curvature of the nanotube [Isj . By con- 
trast to the topography image which is energy integrated, 
dl/dV maps reflect the LDOS at a given energy. It 
clearly displays C-C bond asymmetry, while keeping the 
translational invariance of the lattice. Alternating im- 
age types are observed when going from HOMO-1 to 



LUMO+1 : at -0.78 V, the intensity is maximum above 
bonds of type 1, while at -0.37 V there is a node above 
bond 1 and maxima above bond 2 and 3. At +0.37 
V the pattern is similar to HOMO-1 while at +0.78 V 
it is similar to HOMO. These images demonstrate the 
complementary nature of the wavefunctions at the VHS 
energies, which is completely consistent with the asym- 
metries of the spectra shown in Fig. [IF.. We stress that 
the present measurements were performed far from the 
extremity of the tube and that the topographic image 
presents a hexagonal network with the same height for 
the 3 bond orientations and without any evidence of the 
presence of defects. Such alternating patterns have been 
observed systematically on the semiconducting tubes, ei- 
ther individually lying on the substrate or assembled in 
bundles. 

Our observation are consistent with the general sym- 
metry arguments in the reciprocal space put forward by 
Kane and Mele [lfj. A chemist's approach in real space 
appears to be very fruitful also as we show now. We 
consider first the VHS of a semiconducting tube. The 
corresponding Bloch cigenstates have only a double triv- 
ial degeneracy ±fc, and k lies in the vicinity of the K 
point of the Brillouin zone of the graphene sheet [lj . In- 
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deed, the rolling up of the sheet implies that k should 
be on the so-called cutting lines parallel to the axis of 
the tube and separated by 2tt/L, where L is the length 
of the chiral vector, equal to ird, d being the diameter of 
the tube. If (n, m) are the usual coordinates of the chiral 
vector, the position of the cutting lines with respect to 
the K point depends on the value of n — m modulo 3. In 
the case of a semiconducting tube, n — m = Zp ± 1, and 
the nearest line is at a distance 2ir/3L on the left (resp. 
right) hand side of K when n — m = 3p + 1 (resp. 3p-l); 
see Fig. [21 Close to this point, the energy is proportional 
to \k — K\ and the LUMO state when n — m = 3p + 1 
corresponds to k = K + cjo where go is a vector of length 
2ir/3L pointing along the direction of the chiral vector 
in the negative direction. The other — k value is close to 
the K' point of the Brillouin zone. 




FIG. 2: Brillouin zone of the hexagonal graphene lattice (left) 
with a close-up view at the K point (right). The parallel 
lines are the cutting lines of allowed k for a,n~m — 3p + l 
semiconducting nanotube. The arrow indicates the tube axis 
direction. 

The eigenstates \ip^} are built from two independent 
Bloch states \tp^) = Ci\u{) + C2IU2). These two Bloch 
states are linear combinations of atomic orbitals on 
the two sublattices 1 and 2 of the graphene structure, 
(r\u 1{ 2)) = Er-ei(2) e lk - r " n 4>(r - r^), where 4>(f - r„) is 
the 7r atomic orbital centred on site f* n . It turns out 
that the ratio C1/C2 just introduces a phase factor re- 
lated to the direction of the qo vector, i.e. finally to the 
chiral angle 9, C\jCi — sgn(£7) expi(6 + ip), where cp 
depends on the notations used. For tubes such that 
n — m = 3p — 1, q~o is in the opposite direction and the 
sign of C1/C2 changes. We calculate now the electronic 
density p(r) = \(r\'4'±k)\ 2 ■ This involves an expansion in 
products of overlap functions <j>(f— f^)4>(r — r„). Keep- 
ing the main contributions when — r„ either vanishes 
or is equal to a nearest neighbour vector, we obtain an 
expression for the local density of states n(r,E), which 
is the quantity measured in STS. n(f,E) can be written 
as n{E) x (p a t{r) + Pint(r)), where n{E) is the density 
of states, p a t(r) is the superposition of electronic atomic 
densities and /5int(^) is the interference term responsible 
for the breaking of the sixfold symmetry: 

Pmt(r) - 2<j>(r)4>{r - r Q ) x sgn(E) cos 9 a , 



where 9 a is the angle between the tube axis and the bonds 
Tq : 9 a — 9 or 9± 2ir/3, the chiral angle 9 being here pre- 
cisely defined as the smallest angle between a bond and 
the tube axis, < 9 < tt/6. Because of the dependence 
on the sign of E, the HOMO and LUMO contributions 
are complementary. In the case of the second VHS, the qo 
vector should be replaced by —2qo which shows that the 
singular contribution at the second VHS are also com- 
plementary of those of the first ones. Finally images of 
tubes where n — m = ip — 1 are also complementary of 
those where n — m — ip + 1 since here the qo vector 
shoud be replaced by — qo. Basically two types of images 
are therefore expected at the singularities depending on 
the sign in front of cos9 a in Eq. Q. We call images of 
type I the images corresponding to a positive sign. In 
this case the density above the bonds pointing close to 
the tube axis (angle 9) is reinforced, all the more when 
9 is small, i.e. close to a zig-zag orientation (cos 9 ~ 1). 
The complementary image (image of type II) presents 
stripe reinforcements on the bonds perpendicular to the 
axis (— cos^q, ~ 1/2). Close to an armchair orientation, 
the complementarity concerns principally the bonds close 
to the axis (cos6*q, ~ ±V3/2), the third one being unaf- 
fected (cos0 a ~ 0). As can be seen in Fig. [1] our exper- 
imental results are completely typical of the behaviour 
when cos 9 is close to one. This is consistent with our 
assignement (15,5) for the chiral indices: image of type I 
for the LUMO state when n — m = 3p + 1 [16( . 

This "striping effect" has been simulated in STM im- 
ages for semiconducting nanotubes at bias correspond- 
ing to VHS ± 1 within a tight binding approach [20| . 
To validate further our analysis, we have performed first- 
principles calculations to simulate STM images on semi- 
conducting zigzag (16,0) and (17,0) nanotubes. The DOS 
are reproduced on Fig. [3^,-b. Fig [3j>j show the simu- 
lated conductance images at energies of VHS-2, VHS-1, 
VHS+1 and VHS+2 for successively the (16,0) and (17,0) 
nanotubes. STM simulations reproduce the "striping ef- 
fect" experimentally observed. (16,0) is an — m — 3p+l 
tube and presents a type II behaviour for HOMO and a 
type I for LUMO in total agreement with measurements 
and tight-binding prediction on the 3p+l tubes. Simula- 
tions on (17,0) which isart — m = 3p — 1 tube exhibit a 
complementary behaviour i.e. a type I for HOMO and a 
type II for LUMO. These simulations demonstrate the va- 
lidity of the tight-binding based interpretation presented 
above. 

Metallic tubes We present now STS results on a 
metallic tube. Fig. [4] shows spectra and dl/dV maps as- 
sociated with the first VHS. Remarkably, two kinds of im- 
age types are observed around each VHS. From the tight- 
binding model previously described, the VHS of metallic 
tubes correspond to the four states \ ±K±3qo) instead of 
the doubled degeneracy previously mentioned for semi- 
conducting tubes. An analysis of wavefunction symme- 
tries, similar to the one presented before, do not predict 
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FIG. 3: (a) Total DOS of semiconducting (16,0) and (17,0) tubes, (c) to (f) show a simulated STM conductance images at 
VHS-2 (-0.52 eV), VHS-1 (-0.21 eV), VHS+1 (0.35 eV) and VHS+2 (0.73 eV) for the (16,0) tube, (g) to (j) show a simulated 
STM conductance images at VHS-2 (-0.47 eV), VHS-1 (-0.21 eV), VHS+1 (0.36 eV) and VHS+2 (0.55 eV) for the (17,0) tube. 




FIG. 4: Topographic image (a) of a metallic carbon nanotube (9 — 26 + 1°) recorded during the CITS measurement with 
labelling of the 3 bonds where STS point spectra were measured (b). (c)-(f): Energy-dependent conductance images of a 
metallic tube at -0.81, -0.73, 0.96, and 1.10 V respectively. The arrow indicates the tube axis. 



the two observed complementary wavefunctions. How- 
ever, because of the so-called trigonal warping effect the 
degeneracy between the states corresponding to opposite 
values of 3<fo is lifted (except for armchair nanotubes) 
as observed experimentally and reproduced in calcula- 
tions 21, [22j]. We then expect, as for semiconducting 
nanotubes, complementary images associated with the 
splitted VHS. This is clearly visible in experimental con- 
ductance spectra. Fig. 0fc-f show indeed that bond 1 is 
reinforced at -0.81 V and bond 2 at -0.73 V. Bond 3, 



almost perpendicular to the axis, is less visible in both 
cases, as predicted for a tube close to an armchair orien- 
tation, (9 = 26 ± 1°). At positive bias, on the left hand 
side of the peak at 0.96 V, bond 1 is bright in the dl/dV 
map (Fig. 0k) as expected (image of type I). The conduc- 
tance image shown in Fig. 0f is more complex due to the 
limited experimental resolution, but it exhibits brightest 
spots localized on bond 2 instead of bond 1 in Fig. [4^. 
This is again in line with the tight binding description. 
The observed complementary sequence can be shown 
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to be consistent with the sign of the warping effect: image 
of type I for the lowest VHS at positive energy. To con- 
firm this effect on splitted VHS of metallic tubes, we per- 
formed first principles simulations on a (8,5) nanotube. 
The calculated DOS reported in Fig. [5^ nicely repro- 
duces the splitting of the VHS and also exhibit a dip 
around the Fermi level which will be discussed later. 

The calculated conductance images at the energies of 
the 4 peaks represented in Fig. [5]i-g show alternation of 
images of typel and II confirming the broken sixfold sym- 
metry on metallic tubes. This is totally consistent with 
the observation and tight binding prediction of striping 
effect at split VHS of metallic tubes. 

We conclude with a significant observation of the sym- 
metry breaking at the band edge of the pseudogap. The 
LDOS in (Fig. Eljtop)) measured by STS clearly shows 
a dip around the Fermi level corresponding to the pseu- 
dogap. The dl/dV maps recorded at -30 mV (Fig. [6^) 
and +70 mV (Fig. \Ejp) exhibit complementary symmetry, 
type I at negative bias and type II at positive bias. In- 
deed, this small gap is characterized by a small q perpen- 
dicular to the tube axis. Complementary images of type 
I and II are therefore also predicted in the tight-binding 
model at the lower and upper band edges, respectively, 
in nice agreement with our observations . The first prin- 
ciples of the (8,5) tube also confirms the striping effect 
at the pseudogap as shown in Fig. [Sb-c. The image at 
positive bias is complementary of the image at negative 
bias. Moreover it is complementary of the image calcu- 
lated at the first split VHS peak at positive bias (Fig. 
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FIG. 6: (top) Zoom of the metallic STS around the Fermi 
level, (a) -(b): Energy- dependent conductance images of a 
metallic tube around the Fermi-level at - 30 and + 70 mV. 
The arrow indicates the tube axis. 



[Sf ) confirming the alternation of symmetry observed ex- 
perimentally. 

To summarize, we have shown that STS allows us to 
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image directly the wavefunctions of carbon nanotubes as- 
sociated with band edges and VHS. The wavefunctions 
intensities present patterns with broken sixfold symme- 
try. The images alternate systematically between only 
two types when going from one singularity to another. 
Considering the HOMO orbital, for a 3p — 1 tube the 
electron density is concentrated on the C-C having the 
smaller angle with the tube axis. For a 3p + 1 elec- 
tron density is maximum on the two other bonds. For 
a metallic tube (3p) the sixfold symmetry is expected 
to be recovered. However the trigonal warping lifts the 
degeneracy and pairs of complementary orbitals are re- 
vealed around the energy position of the VHS. These 
results give a complete view of nanotubes molecular or- 
bitals and allow to predict easily any wavefunction of 
any nanotube with a honeycomb lattice. This global pic- 
ture is expected to provide essential insight for all cases 
where the symmetry relation between nanotubes orbitals 
and their surrounding is a key factor such as functional- 
ization, transport devices or nanotubes sorting. 



EXPERIMENTAL AND THEORETICAL 
DETAILS 

Sample preparation and STM setup STM measure- 
ments were performed with a low temperature STM op- 
erating at 5K under UHV conditions (less than 10~ 10 
mbar). Spectroscopy was achieved in the current imag- 
ing tunneling spectroscopy (CITS) mode. Local dl/dV 
spectra were measured with a lock-in amplifier at each 
point of the images. Conductance images as well as local 
point spectra were then extracted from the data. The 
SWNTs were synthesized in a vertical flow aerosol reac- 
tor and deposited in-situ onto commercial gold on glass 
surface previously flashed by butane flame in air. After 
the synthesis, the samples were introduced into the UHV 
system and heated to 150°C during 30 min, for degasing. 
All measurements were performed with tungsten tips. 

DFT and STM simulations Density-functional calcu- 
lations were performed using the SIESTA code within the 
local density approximation as parametrized by Perdew 
and Zunger [231 ]. The core electrons were replaced 



by nonlocal norm-conserving pseudopotentials 24(. A 
double-£ basis set of localized atomic orbitals were used 
for the valence electrons. STM topological images were 
calculated according to the Tersoff-Hamann approxima- 
tion [Hj]. 
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